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ABSTRACT 

A  preliminary  evaluation  of  the  RANS  code  CFX- 
TASCflow  for  predicting  the  hydrodynamic  forces  on 
a  fully  appended  submarine  at  a  Reynolds  number  of 
20  million  is  described.  Except  for  drag,  TASCflow 
makes  good  predictions  of  the  overall  forces  up  to 
moderate  incidence  angles  using  wall  functions  with 
a  371,000  node  grid.  At  high  incidence  angles  the 
predictions  break  down  due  to  inadequate  modelling 
of  separation  and  shed  circulation.  The  predictions 
improve  with  a  finer  grid  and  different  turbulence 
models.  Drag  predictions  improve  when  the  bound¬ 
ary  layer  is  modelled  in  place  of  wall  functions. 


NOMENCLATURE 

Standard  submarine  body  fixed  forces  are  used,  with 
the  origin  located  at  the  center  of  buoyancy. 


K,  M,  N  rolling,  pitching,  and  yawing  moments. 
U  overall  velocity:  \/u2  +  v2  +  w2  . 

u,v,w  body  axis  velocities  of  vehicle. 

X,  y,  Z  axial,  lateral,  and  normal  forces. 
f3  angle  of  drift  (yaw):  tan-1(— v/u) . 

1  INTRODUCTION 

AEA  Technology  and  DREA  are  evaluating  the 
AEA  Technology  Reynolds  Averaged  Navier-Stokes 
(RANS)  code  CEX-TASCflow  as  a  suitable  tool  for 
providing  DREA  with  a  robust  underwater  vehi¬ 
cle  hydrodynamic  force  predictor.  TASCflow  is  a 
second  order,  pressure/velocity  coupled,  finite  ele¬ 
ment  based  control  volume  method  that  uses  a  body- 
fitted,  non-orthogonal,  structured  grid  and  a  multi¬ 
grid  solver.  The  code  runs  in  single  precision  usually 


using  a  k-e  turbulence  model  and  wall  functions,  but 
other  turbulence  models  and  a  two-layer  boundary 
layer  model  are  also  evaluated. 

The  objective  of  this  work  was  to  demonstrate 
how  TASCflow  would  perform  on  underwater  vehicle 
or  model  scale  submarine  problems  using  a  grid  sized 
for  computations  on  a  desktop  workstation  and  built 
by  someone  with  limited  experience.  The  work  was 
conducted  in  three  phases  over  several  years. 

In  the  first  phase,  a  masters  student  working 
with  AEA  Technology  built  a  371,000  node  grid  for 
a  generic  submarine  configuration  consisting  of  an 
axisymmetric  hull  with  a  sail  and  four  symmetrical 
tailfins.  TASCflow  calculations  using  a  k-e  turbu¬ 
lence  model  and  wall  functions  were  performed  for 
the  sub  at  yaw  angles  of  0,  10,  20,  and  30  degrees 
at  a  Reynolds  numbers  of  20  million  (based  on  hull 
length) .  These  parameters  correspond  to  the  test  pa¬ 
rameters  for  experimental  data  that  DREA  and  the 
Institute  for  Aerospace  Research  (IAR)  obtained  in 
the  IAR  9  m  wind  tunnel  in  Ottawa  [lj.  This  allowed 
the  TASCflow  predictions  to  be  validated. 

In  the  second  phase,  wall_  functions  were  again 
used.  First;  the 'grid  nodal  density  in  the  outer  flow 
was  doubled  and  the  calculations  repeated.  Second, 
the  371,000  node  calculations  at  30  degrees  inci¬ 
dence,  where  agreement  with  experiment  was  poor¬ 
est,  were  repeated  with  two  k-u  turbulence  models. 

In  the  third  phiase,  the  phase  1  grid  was  expanded 
into  the  boundary  layer  and  the  phase  1  calculations 
repeated  using  a  two-layer  turbulence  model  and  no 
wall  functions.  The  k-e  model  was  used  in  the  outer 
flow  and. standard  one-equation  and  algebraic  mod¬ 
els  were  used  to  model  k  and  the  turbulence  length 
scale  l  in  the  boundary  layer. 

2  NUMERICS 
2.1  The  Grids 

The  phase  1  structured  grid  is  shown  in  Figure  1.  It 
is  built  from  32  blocks  and  contains  371,000  active 
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nodes.  The  dimensionless  y+  parameter  values,  vis¬ 
cous  sublayer  scaled  measurements  of  the  distance 
of  the  first  layer  of  nodes  off  the  submarine  surface, 
varied  from  1  to  1000  over  the  surface  at  the  var¬ 
ious  incidence  angles.  The  average  y^~  value  was 
163  which  falls  within  the  recommended  range  of 
30  <  y+  <  500  for  the  wall  functions. 

The  grid  had  a  minimum  skew  angle  of  20  degrees 
which  is  the  minimum  angle  recommended  for  the 
TASCflow  code. 

TASCflow  uses  single  precision  floating  point 
arithmetic  to  reduce  memory  requirements  and  acco¬ 
modate  large  problems  on  desktop  computers.  How¬ 
ever,  this  means  that  round-off  error  limits  the  cell 
aspect  ratios  a  grid  can  have.  This  problem  is  al¬ 
leviated  by  the  use  of  wall  functions  which  reduce 
the  need  for  high  grid  densities  near  surfaces  where 
the  highest  aspect  ratio  grid  cells  occur.  The  recom¬ 
mended  maximum  cell  aspect  ratio  for  TASCflow  is 
100:1.  About  6%  of  the  phase  1  grid  aspect  ratios 
were  larger  than  this,  with  the  largest  being  780:1. 

In  phase  2,  the  phase  1  grid  was  built  up  to 
932,000  nodes;  wall  functions  were  again  used.  The 
y+  distribution  was  marginally  improved  and  ranged 
from  6  to  860  with  an  overall  average  of  150.  Aspect 
ratio  and  skew  characteristics  were  unchanged. 

In  phase  3,  the  phase  1  grid  was  expanded  into  the 
boundary  layer  for  a  total  of  835,000  active  nodes. 
The  boundary  layer  was  modelled  using  a  two-layer 
k-e  / k-l  model.  The  grid  had  an  average  of  3  nodes 
within  a  y+  of  15,  a  maximum  aspect  ratio  of  705:1, 
and  a  minimum  skew  angle  of  15  degrees. 

2.2  Boundary  Conditions 

All  boundary  conditions  were  steady  state. 

Either  a  conventional  smooth  wall,  log-law  wall 
function  model  or  the  two-layer  k-e  I  k-l  model  were 
used  on  the  surface  of  the  vehicle. 

Inflow  velocity  (u  =  U cos/3,  v  —  U sin/?,  w  = 
0),  turbulence  intensity  (1%,  to  match  the  wind  tun¬ 
nel),  and  turbulence  length  scale  (about  half  a  hull 
diameter)  were  specified  at  the  two  upstream  outer 
faces  of  the  grid.  When  0  =  0,  this  boundary  con¬ 
dition  was  only  applied  to  the  single  upstream  face. 

The  two  outflow  faces  (one  when  0  =  0)  had 
their  average  static  pressure  set  to  zero. 

An  ‘opening1  boundary  condition  was  used  on  the 
remaining  two  faces  (four  when  0-0)  in  which  the 
static  pressure,  turbulence  intensity,  and  turbulence 
length  scale  were  fixed  at  far  field  values.  The  flow 
direction  was  implicitly  determined  by  the  solution 
and  could  vary  between  inflow  and  outflow  over  a 
face. 
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Figure  2  Convergence  history  of  RMS  residuals  for 
the  371,000  node  grid,  f3  =  20  degrees. 


For  the  phase  1  grid,  a  coalescence  boundary  con¬ 
dition  was  applied  at  the  aft  tip  of  the  hull  where 
the  entire  face  of  the  trailing  grid  block  degenerates 
to  a  point.  This  condition  could  not  be  used  for 
the  large  grids  as  it  prevented  convergence.  Solu¬ 
tions  obtained  without  using  the  coalescence  bound¬ 
ary  condition  had  about  a  10%  variation  in  pressure 
over  this  degenerate  face,  creating  a  minor  disconti¬ 
nuity  in  the  solution  at  the  aft  tip  of  the  hull.  Per¬ 
haps  finer  longitudinal  localized  gridding  here  would 
solve  this  problem. 

2.3  Convergence 

In  most  cases  the  solutions  converged  satisfactorily, 
as  shown  in  Figure  2.  There  were  two  exceptions. 
One  was  the  coalescence  boundary  condition  prob¬ 
lem  discussed  above.  The  other  occurred  for  the  zero 
incidence  calculations  and  was  eventually  traced  to- 
the  multigrid  solver;  with  the  multigrid  solver  turned 
off,  the  zero  incidence  solutions  readily  converged. 

3  RESULTS 

The  results  were  very  good,  as  shown  by  the  body 
axes  force  comparisons  in  Figure  3.*  Except  for  ax¬ 
ial  force  (negative  drag),  which  is  unimportant  for 
vehicle  stability  and  maneuvering  performance  pre¬ 
dictions,  ail  the  forces  were  well  predicted  at  low 
to  moderate  incidence  angles.  This  range  encom¬ 
passes  the  standard  towing  tank  test  range  (0  to  18 
degrees)  over  which  modern  submarines  are  model 
tested  before  being  built. 

Of  increasing  importance  to  submarine  safety 
are  performance  predictions  for  extreme  maneuvers. 

*  Proprietary  considerations  prevent  publication  of  j/-axis 
scales  with  these  data. 


This  requires  knowledge  of  the  hydrodynamic  forces 
at  incidence  angles  larger  than  most  towing  tank  fa¬ 
cilities  can  handle,  such  as  the  30  degree  yaw  an¬ 
gles  obtained  in  the  IAR  9  m  wind  tunnel.  It  is  de¬ 
sirable  to  obtain  this  knowledge  more  economically, 
through  the  use  of  a  CFD  code.  As  Figure  3  shows, 
however,  TASCflow  is  less  accurate  at  predicting  the 
forces  at  high  yaw  angles,  especially  the  important 
out-of-plane  forces  (normal  force  and  pitching  mo¬ 
ment)  which  directly  influence  depth-keeping.  This 
is  attributed  to  the  inability  of  the  code,  grid,  or 
turbulence  model  to  correctly  shed  or  transport  cir¬ 
culation  in  vortices  trailing  from  the  sail  and  hull  at 
large  incidence  angles. 

Figure  4  shows  experimental  total  pressure  and 
transverse  velocity  fields  in  the  propeller  plane  of 
the  submarine,  just  aft  of  the  tailplanes.  The  sail 
tip  trailing  vortex  is  clearly  visible  in  the  experi¬ 
mental  plots  two  hull  diameters  away  from  the  hull 
centerline  at  about  2  o’clock.  The  lower  hull  sepa¬ 
ration  vortex  is  about  one  hull  diameter  away  at  3 
o’clock.  This  lower  vortex  contains  the  image  cir¬ 
culation  from  the  sail  tip  vortex  which,  at  low  inci¬ 
dence  angles,  remains  bound  in  the  hull  afterbody. 
The  out-of-plane  forces  are  generated  when  this  hull- 
bound  circulation  interacts  with  the  crossflow  along 
the  afterbody  (the  Magnus  effect).  These  forces  roll¬ 
off  at  large  yaw  angles  as  the  hull-bound  circulation 
is  shed  in  the  lower  body  separation  vortex. 

Figure  4  shows  that  TASCflow  is  not  modeling 
these  trailing  vortices  well,  which  means  the  trailing 
sail  circulation  remains  near,  if  not  still  bound  in, 
the  hull  afterbody.  Thus,  TASCflow  overpredicts 
the  out-of-plane  forces  at  high  yaw  angles. 

•  ••The  high  incidence  out-of-plane  force  predictions 
improve  when  1)  the  turbulence  model  is  replaced 
with  either  of  two  variations  of  the  k-uj  turbulence 
model  and  2)  the  grid  nodal  density  in  the  outer  flow 
is  increased.  High  incidence  axial  force  predictions 
improve  with  these  changes  as  well,  but  the  two-layer 
model  improves  axial  force  over  the  full  incidence 
range. 

Figure  5  shows  the  TASCflow  total  pressure  con¬ 
tours  and  transverse  velocity  vectors  for  the  large 
grids.  When  the  nodal  density  in  the  outer  flow 
is  increased  (932,000  node  grid)  the  total  pressure 
better  reflects  experiment.  However,  the  transverse 
velocity  vectors  show  that  circulation  is  still  not  ad¬ 
equately  modelled.  When  the  nodal  density  in  the 
boundary  layer  is  increased  (835,000  node  grid)  there 
is  no  improvement  in  the  pressure  predictions. 
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f  371,000  nodes  (wall  functions):  ©  k-e,  a  k-ui,  v  k-u-mod 
x  Exp  data;  TASCnow  j  large  grids  (k-e):  •  wall  functions,  ©  two-layer 


The  371,000  node  calculations  use  wall  functions.  The  large  grid  calculations  (932,000  node  grid  using  wall 
functions  and  835,000  node  grid  using  the  two-layer  model)  use  the  k-e  turbulence  model. 
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Figure  4  Total  pressure  coefficient  contours  and  transverse  velocity  vectors  in  the  submarine  propeller 
plane  at  0  =20  degrees.  The  white  disc  in  the  middle  of  each  figure  is  the  hull.  The  black  annulus  in  the 
experimental  plot  is  a  Bow  region  where  data  was  not  obtained.  The  large  circles  in  the  Bgures  represent  the 
maximum  hull  diameter.  The  velocity  vector  length’s  are  proportional  to  the  transverse  velocity  magnitudes. 


4  DISCUSSION 

These  results  suggest  that  grid  density  in  the  outer 
flow  should  be  increased  and  the  turbulence  model 
optimized  to  improve  modelling  of  shed/ trailing  vor- 
ticity  and  of  the  out-of-plane  forces.  This  is  con¬ 
sistent  with  the  findings  of  other  submarine  CFD 
researchers. 

Zierke  [2]  found  that  about  10  nodes  were  needed 
in  the  core  of  a  vortex  to  prevent  numerical  diffusion 
during  convection.  Figures  4  and  5  lack  this  resolu¬ 
tion  (velocity  vector  density  equals  nodal  density  in 
the  figures),  especially  in  the  sail  tip  vortex. 

The  best  turbulence  model  for  submarine  flows 
has  yet  to  be  determined.  Early  submarine  CFD 
work  [3]  found  that  the  standard  algebraic  Baldwin- 
Lomax  model  required  modifications  to  adequately 
predict  separation  (wall-functions  were  not  used). 


Subsequent  work  [4]  used  the  k-tn  two-equation 
model,  but  it  was  noted  that  k-u>  overpredicts 
Reynolds  stress  in  the  outer  boundary  layer  region. 

AEA  Technology  recently  carried  out  calcula¬ 
tions  with  their  unstructured  CFX-TfC  code  using 
a  334,000  node  version  of  the  phase  1  grid  and  both 
k-e  and  anisotropic  Reynolds  stress  turbulence  mod¬ 
els.  Wall  functions  were  also  used.  The  six-equation 
Reynolds  stress  model  provided  substantially  better 
drag  predictions  than  k-e  but  only  marginally  im¬ 
proved  out-of-plane  force  predictions. 

Future  work  will  address  the  above  issues  in  more 
detail  and  extend  the  validation  to  other  model  con¬ 
figurations  [2].  In  addition,  AEA  Technology  has  re¬ 
cently  provided  DREA  with  an  improved  wall  func¬ 
tion  formulation  and  a  double  precision  version  of 
the  TASCflow  code  for  evaluation. 
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5  CONCLUDING  REMARKS 

TASCflow  does  a  good  job  of  predicting  the  overall 
forces  at  low  to  moderate  incidence  angles,  except 
for  drag.  This  incidence  range  encompasses  the  test  ^ 
regime  of  conventional  towing  tank  facilities.  The 
authors  know  of  no  other  non-experimental  tool  that 
can  match  this  accuracy. 

At  high  incidence  angles,  the  current  TASCflow 
overall  force  predictions  break  down  due  to  inade¬ 
quate  modelling  of  separation  and  shed  circulation. 
These  predictions  improve  with  a  more  optimum  tur¬ 
bulence  model  and  with  finer  gridding  in  the  regions 
of  separated  flow. 

TASCflow  drag  predictions  improve  when  wall 
functions  are  replaced  with  a  boundary  layer  calcula¬ 
tion.  However,  drag  characteristics  are  not  critical  to 
submarine  maneuvering  studies.  In  fact,  the  overall 
results  of  this  work  suggest  that  using  wall  functions 
may  be  an  acceptable  way  to  gain  the  memory  nec¬ 
essary  to  adequately  grid  the  trailing  vortex  fields  in 
the  outer  flow  of  a  maneuvering  vehicle. 
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